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Abstract 

By means of a new multilayer pseudo-spectral moist-convective thermal rotat- 
ing shallow-water (mcTRSW) model in a full sphere, we present a possible 
equatorial adjustment beyond Gill’s mechanism for the genesis and dynam- 
ics of the Madden-Julian oscillation (MJO). According to this theory, an 
eastward-propagating MJO-like structure can be generated in a self-sustained 
and self-propelled manner due to nonlinear relaxation (adjustment) of a 
large-scale positive buoyancy anomaly, depressed anomaly, or a combination of 
these, as soon as this anomaly reaches a critical threshold in the presence of 
moist convection at the Equator. This MJO-like episode possesses a convectively 
coupled “hybrid structure” that consists of a “quasi-equatorial modon” with an 
enhanced vortex pair and a convectively coupled baroclinic Kelvin wave (BKW), 
with greater phase speed than that of dipolar structure on an intraseasonal 
time-scale. Interaction of the BKW, after circumnavigating the entire Equator, 
with a new large-scale buoyancy anomaly may contribute to excitation of a recur- 
rent generation of the next cycle of MJO-like structure. Overall, the generated 
“hybrid structure” captures a few of the crudest features of the MJO, including 
its quadrupolar structure, convective activity, condensation patterns, vorticity 
field, phase speed, and westerly and easterly inflows in the lower and upper 
troposphere. Although moisture-fed convection is a necessary condition for the 
“hybrid structure” to be excited and maintained in the proposed theory in this 
study, it is fundamentally different from moisture-mode theories, because the 
barotropic equatorial modon and BKW also exist in “dry” environments, while 
there are no similar “dry” dynamical basic structures in moisture-mode theo- 
ries. The proposed theory can therefore be a possible mechanism to explain the 
genesis and backbone structure of the MJO and to converge some theories that 


previously seemed divergent. 
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1 | INTRODUCTION 


More than half of the globe is considered “the Tropics”; 
thus a proper understanding of the tropical atmosphere 
is crucial for improving global weather forecasts, as 
well as for projections of global climate change. Mat- 
suno (1966) and Gill (1980) present the foundations of 
equatorial waves by utilizing the one-layer version of the 
shallow-water model on the equatorial #-plane. Let us 
recall that rotating shallow-water (RSW) models are effec- 
tive and low-cost numerical tools with which to study 
the large-scale dynamics of the oceans and atmosphere 
based on vertically averaged primitive equations between 
pairs of material surfaces in vertical coordinates. The 
RSW model in the equatorial f-plane is an archetype 
model for understanding dynamics in the Tropics. The 
derived spectrum of equatorial waves is highly consistent 
with outgoing long-wave radiation (OLR) data (Wheeler 
and Kiladis, 1999), except for a branch corresponding 
to slowly eastward-moving motions associated with the 
Madden-Julian oscillation (MJO). The MJO is the promi- 
nent fluctuation in tropical weather over a 30-90 day 
cycle (Madden and Julian, 1972; 1994). This large-scale 
circulation of enhanced deep convection moves slowly 
eastward from the Indian Ocean and then passes over 
the maritime continent and finally dies out over the 
Pacific Ocean. The MJO is accompanied by extreme rain- 
fall and wind over an area roughly 1,500 km in latitude 
and 4,500 km in longitude (Zhang et al., 2020) and has a 
great impact on both tropical moist-air circulation and 
midlatitude environment variables on the intraseasonal 
time-scale (Vitart and Molteni, 2010). There have been 
several attempts to explain basic features of the MJO in 
terms of various dynamic mechanisms. A comprehensive 
description of the existing state of theories of the MJO has 
been provided in recent reviews by Zhang et al. (2020) and 
Jiang et al. (2020). In particular, Kim and Zhang (2021) 
recently found a slowed-down Kelvin-wave (KW) solution 
that resembles intraseasonal periodicity on the plane- 
tary scale and eastward propagation like that of the MJO 
by adding momentum damping to the linear equatorial 
shallow-water equations in a “dry” environment. Theo- 
ries of the MJO are still diverse and further advancements 
in understanding the crudest dynamics of the MJO are 
unavoidable (Zhang et al., 2020). Previous theoretical 
studies, modeling, and field observations, in terms of 
linear and nonlinear equatorial waves of the MJO, have 
contributed positively to our understanding of the MJO; 
nevertheless, a few dynamical aspects of the MJO are not 
fully understood. Some of these aspects include initiation 
of the MJO, its backbone structure, the role of moist con- 
vection in determining the phase speed, the evolution of 
potential temperature in the lower and upper troposphere 


during MJO evolution, and the roles of Rossby and Kelvin 
waves in the MJO’s eastward propagation mechanism. 
The present study aims to shed light on these aspects 
and to show that there are high degrees of convergence 
between some theories that previously seemed divergent. 
The theory that is particularly relevant to this research is 
one that regards the MJO as a “hybrid structure”, which 
consists of a “quasi-equatorial modon” and convectively 
coupled baroclinic Kelvin wave (CCBKW: Rostami and 
Zeitlin, 2020a). This theory has been developed in a step- 
wise manner, with a hierarchy of models with different lev- 
els of complexity, in such a way that each step has provided 
a solid basis for further theoretical studies in the next step. 
All steps are described below for better understanding. 


1. In the first step, it was shown that the barotropic equa- 
torial modon has a persistent eastward propagation for 
a long time span in a “dry” or adiabatic environment 
(Yano and Tribbia, 2017; Rostami and Zeitlin, 2019a). 
A modon is defined as an exact, form-preserving, uni- 
formly translating, horizontally localized, nonlinear 
solution to the inviscid quasigeostrophic equations. 
Yano and Tribbia (2017) illustrated a barotropic modon 
from solutions of Tribbia (1984) and Verkley (1984) 
on a sphere. Rostami and Zeitlin (2019a) obtained a 
similar coherent dipolar structure, the so-called equa- 
torial modon, via the process of ageostrophic adjust- 
ment from the analytical asymptotic modon solu- 
tions of the vorticity equation obtained at the limit of 
small pressure perturbations in the equatorial f-plane 
(see the Appendix). Numerical simulation with a 
high-resolution shallow-water model in the equatorial 
B-plane that is initialized with analytical asymptotic 
modon solutions leads to a steady dipolar structure, 
which propagates eastward coherently. One fundamen- 
tal difference of our proposed MJO theory encompass- 
ing the modon as one of the key components of its 
proposed structure is that equatorial modons have an 
exact asymptotic solution and crystal-clear structure; 
furthermore, eastward propagation and small phase 
speed are innate characters of them that exist even in 
“dry” environments. Such a feature does not exist in 
moisture-mode theories, in which no dry dynamical 
basic structure exists. The applied dynamical regime 
in the RSW model by Rostami and Zeitlin (2019a) cor- 
responds to the so-called long-wave approximation in 
oceanography (Gill, 1982), which arises at the limit of 
small pressure variations and gives, to leading order, 
a vorticity equation that does allow for modon solu- 
tions. This quasibarotropic regime is the RSW analog 
of the Charney (1963) nondivergent balance model 
for the equatorial atmosphere. According to scale 
analysis by Charney, synoptic-scale tropical adiabatic 
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motions should obey barotropic vorticity dynamics. 
The modon solutions in a quasigeostrophic (QG) model 
with a finite deformation radius on a midlatitude tan- 
gent plane were already introduced by Larichev and 
Reznik (1976) and then Flierl et al. (1980) showed the 
effects of baroclinicity with their two-layer generalized 
version of the QG model. Nevertheless, the QG approx- 
imation loses its full validity in the vicinity of the equa- 
torial £-plane. It is therefore no wonder this approach 
has been overlooked by the community. In the first step, 
the main downsides of the equatorial modon’s struc- 
ture compared with the MJO were incompatibility in 
vertical structure due to the barotropic structure of the 
equatorial model and lack of a mechanism to explain 
its genesis. The presence of dipolar structure in the 
MJO is well established in the literature (e.g. Lau and 
Phillips, 1986; Matthews, 2000; Zhang and Ling, 2012); 
indeed, Wang et al. (2019) analyzed the vorticity vari- 
ability associated with the MJO and reported the exis- 
tence of a dipolar-vortex signal. This analysis supported 
the existence of a strongly nonlinear modon solution. 

. In the second step, it was demonstrated, for the 
first time, that the equatorial modon’s structure can 
also emerge from the process of equatorial adjust- 
ment of a localized large-scale depression-type dis- 
turbance as soon as it reaches a triggering state in 
the standard moist-convective rotating shallow-water 
(mcRSW) model on the equatorial f-plane (Rostami 
and Zeitlin, 2019b). The critical threshold, triggering 
state, or tipping point in this context refers to a state 
in which a small additional perturbation transmits the 
large-scale evolution and structure of the system to an 
eastward propagation. From this step on, there is a tran- 
sition from analytical to numerical solutions. Notably, 
these slowly eastward-propagating modons arise from 
such equatorial adjustment only in the presence of 
moist-convective feedback, not in the “dry” case. How- 
ever, once an equatorial modon is created, it can survive 
for a long time span even in a “dry” environment. This 
finding opened a new window to look for the dynam- 
ical relation between equatorial adjustment and stim- 
ulated eastward-propagating structures. Indeed, this 
result undermined the universality of the well-known 
classical Gill (Gill, 1980; 1982) theory of tropical cir- 
culation due to localized heating. According to Gill’s 
mechanism, when the imposed heating is symmetric 
about the Equator, the response can be interpreted as 
steady, damped Kelvin and Rossby waves extending to 
the east and west of the heating, respectively. Zhao 
et al. (2021) could produce a similar eastward propaga- 
tion due to geostrophic adjustment like that of Rostami 
and Zeitlin (2019b) with a moist-convective RSW model 
on a sphere. 
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3. In the third step, further studies showed that increas- 


ing baroclinicity eventually leads to loss of coherency 
and arrest of eastward propagation (Rostami and 
Zeitlin, 2020b). Following the process of equato- 
rial adjustment, Rostami and Zeitlin (2020a) studied 
the process of relaxation (adjustment) of localized 
large-scale pressure anomalies in the lower equato- 
rial troposphere, and showed that such phenomena 
can engender coherent structures strongly resem- 
bling some of the crudest features of MJO events, as 
seen in vorticity, pressure, phase speed, time span, 
and condensation patterns. It is worth noting that 
by using a Lagrangian atmospheric model, Liang 
et al. (2021) re-examined the hypothesis of the gen- 
eration of MJO structure due to equatorial adjust- 
ment of large-scale localized heating (Rostami and 
Zeitlin, 2020a) by imposing warm sea-surface temper- 
ature anomalies of different aspect ratios and mag- 
nitudes and could observe robust generation of an 
MJO event. The aforementioned equatorial adjustment 
simulation in moist-convective environments by the 
mcRSW model was accomplished by forcing a nega- 
tive pressure anomaly, with an invariant gradient of 
potential temperature. 


. In the present study, we intend to advance the recently 


proposed theory by Rostami and Zeitlin (2020a) one 
step further by completing the missing part of their 
model, which is the variation of potential temperature. 
In order to overcome this bottleneck, we have imple- 
mented a multilayer pseudo-spectral moist-convective 
thermal rotating shallow-water (mcTRSW) model (see 
Section 2) to capture evolution of the buoyancy field 
and its interaction with other fields explicitly in 
the lower troposphere. All the principal equations 
of the mcTRSW model and numerical scheme in 
this study are different from those of the mcRSW 
model by Rostami and Zeitlin (2020a). Our results 
show that an MJO-like skeleton can be agitated in a 
self-sustained and self-propelled manner as a “hybrid 
structure” right after a large-scale localized heating 
and/or depression-type disturbance in the lower tropo- 
sphere reaches a triggering state in a moist-convective 
environment. The “hybrid structure” is constituted by a 
combination of a dipolar structure or quasi-equatorial 
modon and a CCBKW that lasts for an intrasea- 
sonal time-scale. The baroclinic Kelvin wave (BKW) 
then detaches from the dipole due to its higher east- 
ward phase speed. It is shown that catching up and 
interaction of the circumnavigated BKW with a new 
large-scale positive buoyancy anomaly may trigger the 
recurrent generation of a new MJO-like structure, and 
its interaction with the residue of the eastward-evolving 
part from the previous cycle weakens the residue of 
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the eastward-propagating wave. The latter interaction 
is accompanied by subsequent destabilization due to 
shear in vorticity and an elongated temperature front. 
The proposed mechanism of this theory captures a few 
basic features of the MJO structure as seen in its east- 
ward propagation with low phase speed, quadrupolar 
structure, convective activity, condensation patterns, 
vortex-pair structure, westerly and easterly inflows in 
the lower and upper troposphere, and episodic recur- 
rent generation on an intraseasonal time-scale in the 
presence of a large-scale warm pool. 


The article is organized as follows. In Section 2 we 
describe the multilayer pseudo-sparse spectral mcTRSW 
model with minimal parameterization in the full sphere. 
Then, in Section 3, by means of the mcTRSW model we 
examine the adjustment scenarios by applying the model 
for the equatorial adjustment of large-scale localized posi- 
tive potential temperature and negative thickness anoma- 
lies. We show that both approaches coincide in generat- 
ing the “hybrid structure” with a high order of similarity. 
Section 4 contains a comprehensive comparison of the 
results with the MJO. Conclusion and further discussion 
are in Section 5. 


2 | MODEL DESCRIPTION AND 
SETUP OF NUMERICAL 
EXPERIMENTS 


2.1 | Multilayer thermal rotating 
shallow-water model 


In spite of remarkable theoretical frameworks to affirm 
the capabilities of the TRSW model, a sporadic number of 
multilayer numerical implementations of mcTRSW mod- 
els on a sphere exist among atmospheric models. In the 
classical setup of shallow-water models, fluid is homoge- 
neous, incompressible, and under hydrostatic balance (cf. 
Zeitlin, 2018), so gradients of mean temperature and den- 
sity are invariant in this setup. However, the TRSW models 
employ inhomogeneous layers, which allow horizontal 
variations in material properties in the RSW model (Schopf 
and Cane, 1983; Ripa, 1993). Thus, the TRSW equations 
are also called the inhomogeneous-layer model in the liter- 
ature (cf. Ripa, 1996). The concept of TRSW has a profound 
theoretical background and a relatively long history in the 
literature. For instance, this approach has been applied 
to examine dynamical effects of a well-mixed planetary 
boundary layer (Lavoie, 1972) or to formulate a reduced 
two-dimensional system governing deep atmospheric 
motions under simple classes of stratification (Salby, 1989). 
The TRSW model in some cases has been explained as 


a model for an upper active layer of fluid on top of a 
lower inert layer (Warneford and Dellar, 2013). Ripa (1993) 
reorganized a framework for a multilayer thermal RSW 
model and a low-frequency approximation (Ripa, 1996) 
of the model. Recently, another generalization of Ripa’s 
single-layer model (Ripa, 1995) to an arbitrary number 
of layers with stratification and shear has been rederived 
by Beron-Vera (2021). Using a novel well-balanced cen- 
tral upwind scheme, Kurganov et al. (2021b) compared 
the TRSW model with the isothermal RSW model for 
dynamical processes of evolution of isolated vortices in 
the midlatitude #-plane. In another approach, Zerroukat 
and Allen (2015) presented a TRSW model that is derived 
from a three-dimensional Boussinesq approximation of 
the hydrostatic Euler equations, where they used the abso- 
lute temperature as the thermodynamic variable. 

In the first place, we include variation of potential tem- 
perature in the “dry” multilayer TRSW model, in which 
the density of a fluid can vary. For this purpose, we 
start with vertical integration of the atmospheric primitive 
equations in terms of the pseudo-height isobaric verti- 
cal coordinates (Hoskins and Bretherton, 1972). We write 
them down in Cartesian coordinates (x, y, Z): 


m + V - Vvh + fÈ X vp = -Vh Ë, (1a) 
oP 0 

— =g— 1b 
az 80," (1b) 
vort 6 (1c) 

OZ 

00 00 

a , —=0. 1d 
ay i Val ew (1d) 


These equations present the horizontal momentum 
(Equation 1a), the hydrostatic relation (Equation 1b), the 
incompressibility (Equation 1c), and the conservation of 
potential temperature (Equation 1d) respectively. Here 
VY = (u,v, w) stands for the 3D velocity field, py = (u, v) is 
its horizontal part, ® is geopotential, V and V, are three- 
and two-dimensional gradient operators, respectively, 0 is 
the potential temperature of the fluid, with a reference 
potential temperature 0, at sea level, and g and Z stand 
for gravity acceleration and the unit vector in the vertical 
direction, respectively. By considering an m-layer configu- 
ration, we proceed with the standard procedure of vertical 
averaging of Equation 1a between the pseudo-height sur- 
faces z(x,y, t) and z11(x,y,t) and make the mean-field 
(columnar motion) approximation, which allows one to 
split averages of products into products of averages. We will 
count the layers from bottom to top in this section, that is, 
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Zo < Z1 < +++ < Zm, hi = Zi — Zi-1, Where m is the number 
of layers and zo is the lowermost material surface. We thus 
obtain the following “master” equation: 


(Zi — Zi-1) (> + (vi : Vwi + fÈ x vi) 


ot 
1 26; 
=-V Gx —Z4) + aeon ) 
2 b; 
+ (zi) Vzi — B(Zi-1) VZi-1. (2) 


After averaging, we omitted the subscript “h”. 
Velocity v and V will be understood from now on as 
two-dimensional vectors. We will use the Lagrangian 
derivative d;/dt = 0/ðt+v;i- V. The layer-averaged 
buoyancy variable in the atmospheric case is defined 
as bi = 20;/ Os, where 6; is the layer-averaged potential 
temperature. 

From the hydrostatic equation (Equation 1b) we get 


D(z) = O(Z-1) + hibi = P(Zo) 


bo i 
89; 
+ DI hy = Po) + X hyb, (3) 
j=l “S j=l 


where j is the index of the sigma notation. For each layer of 
the TRSW model, the corresponding vertically integrated 
equations can be written as 


divi x 1 
ae +f? xvi = b;iVzi + 5 hiV by 
i 
-V (ow + Yn, j (4a) 
j=l 

oh; 

— + V(hjv;) = 0, (4b) 
ot 
dibi 
— =0. 4 
T (4c) 


Potential temperature should increase with (pseudo-) 
height in the atmosphere: namely, the layer with the high- 
est b is at the top, which is upside down relative to the 
oceanic model. The pseudo-height is a function of pres- 
sure: that is, it is a modified pressure coordinate, so it can 
vary on the ground (the lowermost material surface, Zo) 
as a free surface while the geopotential ®|,, is constant. 
A “flat top” boundary condition with constant pressure 
and, as a consequence, constant z is usually imposed at 
the uppermost material surface, which can be roughly 
interpreted as the tropopause. In order to get a two-layer 
TRSW reduction of the atmospheric primitive equations, 
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we introduce the thicknesses of the layers: hı = Z1 — Zo, 
hz = Z2 — Z. The boundary conditions are ®|,, =constant, 
Vz = 0. After replacing the pseudo-height surfaces in 
Equation 4a with thicknesses and taking into account the 
boundary conditions, the momentum equations in final 
form, after some simplifications, become 


a +f? xv = —ShiVbi — bi V(hy + hə), (Sa) 
dzv2 A 1 
Ft Lex va = 5MVb — V (nbs + hab). (5b) 


2.2 | Thermo-quasigeostrophic balance 
relations for the barotropic and baroclinic 
state 


To obtain the thermo-quasigeostrophic (TQG) limit of the 
“dry” TRSW model in the barotropic case without topog- 
raphy, we utilize the following scaling: 


(u,v) ~ U, (x,y) ~ L, h ~ Ho(1 + An), 
b ~ B(1 + 24b"), f~ PA + Ay), (6) 


where A is a small nondimensional parameter and 7 is a 
deviation from thickness at rest Hp. The nondimensional 
momentum equation, after ignoring the higher-order 
terms on the right-hand side, is 


dv 


Ro 
dt 


++ yexv= Syne), (7) 
where Bu = BH /f?L? is the Burger number. With respect 
to the standard for the quasigeostrophic regime, we 
assume that deviations of thickness and buoyancy from 
their mean values are small, e.g. the same order as the 
Rossby number Ro = U/foL = A. If we assume Bu ~ 1, 
then Equation 7 becomes 


Z2xv=-V(n+b’) = Vy, (8) 


where y = 7 + b' is the geostrophic streamfunction. 

By similar scaling, we can write down the correspond- 
ing nondimensional momentum equations of a two-layer 
TRSW model that can be deduced from Equation 5. We 
consider a constant geopotential at the lower boundary. We 
use the following scaling and represent thicknesses and 
buoyancy fields as 


x,y) i L, 
hy2 ~ Ho(ô12 + Amz), 
bi2 ~ Bo(Bi2 + 4b} ,), 


(1,2, V1,2) ~ U, 


f ~ fo; (9) 
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where 61.2 = Hi2/Hpo denotes the nondimensional thick- 
ness, 71,2 the deviation of thicknesses from 6}, biz the 
buoyancy perturbations, and indices 1,2 represent the 
lower and upper layer respectively. This scaling yields 


Ro + sin(@)Z x vz 


1 
Ieo (Sve, + By Vim + n) ; (10a) 


Ro? + singe x vz 
dt 
ô 
ze (avs + + Vb, +BiVm +B.¥n ) , (10b) 


where ¢ is the latitude, so the thermo-geostrophic balance 
relations between velocity and gradients of pressure and 
buoyancy perturbations are 


sin(p)Z x vı = “Vb, + By V(m +m), (11a) 


sin(p)Z X v2 = 6, Vb; + vy, + BıVnı + B72. (11b) 


2.3 | The multilayer mcTRSW 

Let us recall that phase transitions of water vapor with the 
related latent heat release have already been included in 
the so-called mcRSW model (Bouchut et al., 2009; Lam- 
baerts etal., 2011a; 2011b; 2012). Evolution ofthe vertically 
averaged bulk of humidity and the corresponding effects 
of condensation and latent heat release upon the air col- 
umn are defined in the mcRSW model in a simple albeit 
consistent way. The mcRSW model follows up the semi- 
nal ideas of Gill (1982). An improved version of the model 
can include precipitable water, vaporization, and precip- 
itation (Rostami and Zeitlin, 2018). The mcRSW model 
has already been applied to study the effects of moist 
convection on dynamics of large-scale Earth and plane- 
tary jets and vortices (Lambaerts et al., 2012; Lahaye and 
Zeitlin, 2016; Rostami and Zeitlin, 2017; 2020c; 2022; Ros- 
tami et al., 2017). Recently, Kurganov et al. (2020) showed 
that the one-layer configuration of the mcTRSW model, 
with a well-balanced central upwind scheme, could cap- 
ture horizontal gradients of potential temperature and 
their evolution due to moist convection in a self-consistent 
way. Other well-balanced schemes have already been 
applied for thermal shallow-water equations (Chertock 
et al., 2014; Sanchez-Linares et al., 2016). 

Inclusion of moist convection in this study is an exten- 
sion of the one-layer mcTRSW (Kurganov et al., 2020; 
2021a). We start by integrating the Lagrangian conser- 
vation of the linearized equivalent potential temperature 


(moist enthalpy) in the primitive equation of a pair of 
pseudo-height material surfaces, Z;-ı and z;: 


ee ed £ 
lim E (o T Ea) = s| dz, (12) 


Zi +e 
where d/dt(...) = ð...) + V-V(...) is the 
three-dimensional (3D) Lagrangian derivative, V = 
(u,v,w) stands for the 3D velocity field, 0; represents 
potential temperature, qi represents specific humidity, 
£ is the specific heat of vaporization, Cp denotes the 
specific heat at constant pressure, and S; is an external 
thermal forcing. By considering vertical velocities of the 
fluid as w; = dz;/dt + Wj, wi-1 = dz_1/dt + W;-1, where 
W; presents the extra vertical velocity across the material 
surface i due to convective fluxes, and the columnar bulk 
of humidity as lim,_,o L aed Cp) qi = qi, we obtain 


AGin) + AG + (24% )oy-Way=5. 03 


i 


where Aj(...) = 0/dt(...) +V-vi(...), limeso et = 
Si, vi is the vertically averaged horizontal velocity in 
layer i, and 0; stands for the vertically averaged hori- 
zontal potential temperature. Indeed, we ascribe the val- 
ues of dependent variables 6(z;), 0(Zi-1), qZ), q(Zi-1) at 
the vertical boundaries to their average values inside the 
layer. We therefore obtain lim,_.o[0(z — €) = 0(z-1 + €)] = 
6; and similarly lim, o[(£/Cp)q(zi — €) = (L/Cp)q(Zi-1 + 
€)] = qi/hi. Boundary values can vary with different con- 
figurations. In this study, we are looking for the effect of 
dominant nonadiabatic processes and buoyancy with min- 
imal parameters of convective processes. Thus we consider 
just the latent heat release due to condensation and do 
not include precipitatable water and related vaporization 
in the derivation of a simplified two-layer system with a 
rigid lid. Note that there is no prognostic distribution of 
water vapor in the model and that its initial distribution 
is uniform. Condensation, C;, leads to heating and conse- 
quently enhances the total, vertically integrated, potential 
temperature in the corresponding layer. The potential tem- 
perature is proportional to the specific entropy and hence 
to the heating, with some proportionality to constants 
y,y* > 0, which gives 


A(6:h;) = Ci + "S; (14) 


where C; denotes the latent heat release due to conden- 
sation, which in this simplified model occurs just in the 
lower layer. Therefore, y as the coefficient of C; shows the 
contribution of C; due to its latent heat release as a source 
for the bulk (vertically integrated) potential temperature 
and hence the bulk entropy of the air column. Similarly, 
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y* indicates the contribution of S; to the bulk entropy 
of the air column. y and y* can vary from layer to layer. 
Indeed, we apply S; = R; as a simplified thermal Newto- 
nian cooling, R; = —(hjb; — HiBi)/tr, where r, is relaxation 
time. Substituting Equation 14 in Equation 13, rescaling 
qi = (£g/Cps) f qdz and Ri = (g/0s)Ri, and dropping the 
bars from now on results in the following equation for a 
two-layer system: 


1 


Ahe I 
by +q/h 


[a-n-O+@*-1) Si], a5) 


where C = C, is condensation in the lower layer. From 
Equations 15 and 14, we obtain the Lagrangian derivative 
of the buoyancy in the lower layer: 


dibr =0,b, +v: Vb; 


dt 
aa (ees =) 1) (-C) 


hy bı + (qi/hi) 
+ 1) sı , (16) 


2 (qer -1) 
bı + (qi/hi) 

Consequently, for the upper layer, by considering 

Aq2 = C and A(bzh2) = y* S2, we obtain 
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The magnitude of q;/h; with respect to other terms like 
b; is not negligible, so to capture the precise effect of water 
vapor it has been kept in the model. One interesting feature 
of Equations 15 and 16 is that, ify* = 1, then S; will disap- 
pear on the right-hand side of the h, equation and appear 
just on the right-hand side of the buoyancy equation like 
the traditional form. Figure 1 represents the sensitivity of 
the model to y* with respect to the energy density of each 
layer, which is defined as 


B= | E (hw + bh?) dx dy. (19) 


Greater values of y* yield less loss (gain) of energy density 
in the lower (upper) layer. 

The convective component of vertical velocity pro- 
duces Stokes drag, which is included in the model 
(Equation 20a) similarly to Lambaerts et al. (2011b). 
One interesting aspect of the mcTRSW model is that any 
external heating (S;) affects both potential temperature 
and layer thickness, as in atmosphere thermodynam- 
ics. For instance, by applying Newtonian cooling, we 
see this two-way effect in the following Equations 20. In 
comparison, the heating effect in typical parameteriza- 
tion in modeling the Earth or Mars atmosphere (Rong 


Ah: = —1 [ (-C)+(1-7) S| (17) and Waugh, 2004; Warneford and Dellar, 2013; Seyon 
bz + q2/ħ2 et al., 2017; Rostami et al., 2018) appears either in the 
equations of conservation of potential temperature or 
dob in the incompressibility equation, but not both of them. 
dt =rba + Va» Vb2 Therefore, all-inclusive, the two layer mcTRSW model 
1 q2 (q2/h2)(y* — 1) reduction of the atmospheric primitive equations with 
The. l- hy (=C) + bz + (q2/h2) +1) Sa). minimal adjustable parameters, assuming b; >> qi/hi, 
(18) if any, and inclusion of bottom topography (hp), can be 
0.1 yal 0.1 Tar "a 0.1 xed 
= E- Pe ae nae 
i” ce) Perera eo 
o 0.05 0.05 a 0.05 a 
3 = go aP Jerre aede 
ù 0 ù 0 H own 
È 2 2 
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FIGURE 1 Evolution of the thermal shallow-water energy density, as seen in its deviation from the initial value, in the lower layer (left 


panel), upper layer (middle panel), and overall (right panel) during equatorial adjustment due to an initial positive buoyancy anomaly. Ep is 


the initial value of the energy density [Colour figure can be viewed at wileyonlinelibrary.com] 
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defined by the following set of equations: 


Te +f2 xv = -h Vbi — bıV (h» +hı +h), (20a) 
dv r 1 
aE + fz X V2 = -32Vb2 =, V (hıbı) > b2V(hp + h2) 
l=y 
_ = C 20b 
a (v2 — vı) (C), (20b) 


hy — (By, /b1)Hy 


Tr 


1 = 
dihi +V - (hwi) = a Cora-/) 
(20c) 


hy — (B2/b2)H 
dihz + V - (hava) = B ao +a -pp E/H 


r 


(20d) 
bı — (H,/h,)B 
ihare The ug- BO (20e) 
hı Tr 
b> — (H> /h2)B 
TEE R E E A (208) 
hy Tr 
0914+ V -(qivi1) = Es — C, (20g) 
0:92 + V - (qvı) = C. (20h) 


The relaxational parameterization of condensation, C, 
with a relaxation time, 7, is 


C 


ga - Zu- 21) 


where H denotes the Heaviside (step) function and E; the 
bulk formula for surface evaporation: 


|v] 


[Vmax| 


Es =a 


(Q? — qi H(Q — q1), (22) 


where |Vmax| denotes the maximum amplitude of the 
lower-layer velocity field, QS is the saturation thresh- 
old of bulk humidity, and «œ is an adjustable coeffi- 
cient. In the simplified version, the saturation value QS 
is taken to be constant, but it can depend on poten- 
tial temperature according to the Clausius—Clapeyron law 
(cf. Yano et al., 1995). Parametrization of the conden- 
sation (Equation 21) is of the Betts—Miller type (Betts 
and Miller, 1986), and the parametrization of evaporation 
(Equation 22) is standard in modeling ocean-atmosphere 
exchanges (Katsaros, 2001), except for renormalization by 
the maximum velocity, which is included for convenience. 
To respect the minimal parametrization approach, we have 
included just radiative relaxation and surface evaporation 


(Equation 22) as the only influences of the boundary layer 
in the model. An explicit minimal Newtonian viscosity was 
added to the equations of momentum and buoyancy of 
the model to inhibit the development of small-scale con- 
vective numerical instabilities. Derivation of the spectrum 
of linear waves in the mcRSW model in the long-wave 
approximation is discussed in Rostami and Zeitlin (2020a). 

It is worth mentioning that a wind-induced surface 
heat exchange (WISHE) moisture mode in a simplified 
analytical model produces eastward-propagating modes, 
with the growth rate greatest at the largest wavelengths 
(Fuchs and Raymond, 2005; 2017; Fuchs-Stone, 2020). 
In the convective parameterization of WISHE, variation 
of moist entropy over time is directly linked to zonal 
velocity via the WISHE parameter, which is positive for 
imposed mean westerlies and negative for mean easterlies 
(Fuchs-Stone and Emanuel, 2022). Therefore, perturba- 
tions in surface wind speed enhance the surface enthalpy 
flux, which consequently increase the moist convective 
mass flux and temperature. 


2.4 | Pseudo-spectral numerical method 
and setup 


Numerical methods for the experiments are devised 
according to the Dedalus algorithm, which employs 
spin-weighted spherical harmonics (Lecoanet et al., 2019; 
Vasil et al., 2019). Gelfand and Shapiro (1956) intro- 
duced spin-weighted spherical harmonics for the first 
time in their Lorenz group studies. One advantage of 
spin-weighted spherical harmonics is that a combination 
of spin-weighted spherical harmonics and spinor basis 
vectors renders differentiation on the sphere similar to a 
Fourier series; that is, diagonal wavenumber multiplica- 
tion that remains regular everywhere. We can therefore 
dump traditional singular gradients at the poles. Spectral 
discretizations are specified by the Fourier spectral basis 
for each spatial dimension. It is worth mentioning that 
some other methods have also been proposed to over- 
come problems with coordinate singularities in spheres 
and discs (Boyd and Yu, 2011). The schematic picture of 
a spin-weighted spherical harmonic transform for a set of 
shallow-water equations is 


0u+gVh+2Qe,xu=F-u-Vu, (23a) 
dh + HV -u = F(h) — V - (hw), (23b) 
where u is a 2D velocity field in spherical coordinates, 
U = Upep + Upeo, Where @9,e, are orthonormal coor- 


dinate unit vectors on the surface of the sphere, and 
0<0<xr and 0< g< 2r represent polar and circular 
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angles respectively. Velocity fields and unit vectors on a 
spin-weighted basis are 


1 l 1 l 
U, = —(uọ tity), €, = ie + leg). (24) 


2 2 


We therefore obtain the transformed equations of 
Equation 23 as 


OU, + 8V4h + i2Q cos(0) u, = F, — [u.V uly, (25a) 
du- + gV_h — i2Q cos(0) u- = F_ — [u.V u]-, (25b) 


ðh + H(V u- + V_u,) = Fh) — V (hu) — V_(hu),, 
(25c) 


where V, = (/V2(Vo +i Vọ). For a full description of 
the calculus for vector and tensor operations on the unit 
2-sphere with spin-weighted spherical harmonics and its 
comparison with traditional spherical harmonics, see Vasil 
et al. (2019). By specifying a dealiasing scale factor of 
3/2, we set a 384 by 768 global grid along the latitudi- 
nal and longitudinal directions, respectively. The two-layer 
mcTRSW is calibrated with barotropic equatorial Rossby 
deformation radius Lg = (sH /B)} 2, where H = H; + H, 
is the total nonperturbed layer depth, Hı = (1/3)H, H2 = 
(2/3)H in this study, and p denotes the gradient of the 
Coriolis force in the meridional direction. The zonal and 
meridional velocity, time, and Earth radius scales are as 
follows: 
L~ La, 


(u,v) ~ pE, t~ a a~ La. (26) 


d 


The adjustment experiments were initialized with a neg- 
ative thickness anomaly or positive buoyancy anomaly 
in the lower layer with a/La = 1.5. Other ratios of 


—Time =0 
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a/La = [1.0, 1.50, 1.72,1.88] have been tested as well to 
examine the robustness of the results. The initial nondi- 
mensional thickness anomaly is, according to a simplified 
a-Gaussian equation, 


a 1/0 o 
kW =h-H=h V2e=—r(=+-) G(£,-45), 
(o o 2 Oo 


where r is the spherical distance r= D(Lato, Long, 
Lat, Lon) = D(0, 0, A/Tras, ®)/Ln, Fas determines the aspect 
ratio and Ly affects the spatial perturbation sizes. We set 
Ly = 0.25, ras = 1.8,0 = 4. The method of calculating r is 
based on great-circle distance with respect to the reference 
point according to the Haversine formula. r(x) denotes 
the Gamma function, G(r, a) = [1/T"(a)] f, “e't*—ldt, while 
h and o determine the amplitude and the steepness of 
the thickness anomaly, respectively. We set the maximum 
amplitude of the perturbation as —0.07H for experiments 
that were initialized with a negative thickness anomaly 
in the lower layer and a similar disturbance with positive 
sign for the buoyancy anomaly. Water vapor is initial- 
ized uniformly and close to the saturation value, Q5, in 
the lower layer, Qo = 0.99Q$, and far from the saturation 
value in the upper layer. We set y = 0.15 and y* = 2 for 
the experiments. 

The model produces converged solutions to higher 
precision by applying higher resolution. Propagation of 
eastward fronts of buoyancy and condensation fields 
due to events of a large enough buoyancy anomaly in 
the lower layer, (b,), shows how the model overcomes 
(sub)mesoscale instabilities due to thermal effects in the 
presence of nonlinearity of moist convection (Figure 2) 
and captures such extreme events. While the package of 
buoyancy anomaly pursues the condensed air zone, the 
frontal gradient of the buoyancy anomaly increases with 
enhanced convection until it reaches a ceased state of east- 
ward propagation and damped condensation. We tested 


—Time =0 


-50 0 50 100 150 
Lon 


Snapshots of buoyancy and condensation in the lower layer at 0.2° latitude during equatorial adjustment of an initial 


positive buoyancy anomaly at time = 0, 14, 29, 43, 58 [day], consecutively from left to right (y* = 1) [Colour figure can be viewed at 


wileyonlinelibrary.com] 
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Contour levels of baroclinic thickness anomaly (colors) and corresponding spherical velocity field (arrows) at initial stages 


of equatorial adjustment due to an initial negative thickness anomaly in “dry” (upper row) and moist-convective (lower row) environments. 


Time = 1,4, 13, 24 [day], consecutively from left to right [Colour figure can be viewed at wileyonlinelibrary.com] 


the model for the benchmark problem of Gill’s mechanism 
due to a depressed thickness anomaly, as presented in the 
upper panel of Figure 3. The results are fully consistent 
with those of the RSW model with finite-volume scheme 
(cf. Rostami and Zeitlin, 2020b). 


3 | RESULTS OF EQUATORIAL 
ADJUSTMENT 


3.1 | Equatorial adjustment beyond 
Gill’s heat-induced tropical circulation 


Recently, it has been found that the well-known Gill 
theory of tropical circulation due to localized heating is 
not universal (Rostami and Zeitlin, 2019b). According to 
Gill’s theory, large-scale localized heating on the equato- 
rial f-plane leads to the generation of westward-moving 
Rossby and eastward-moving KWs, while Rostami and 
Zeitlin (2019b) illustrate two other scenarios exist that are 
missed in Gill’s theory. Both scenarios occur right after 
the pressure anomaly or potential temperature anomaly 
reaches a tipping point in the presence of efficient nona- 
diabatic processes like moist convection. These scenarios 
are (i) vortices traveling respectively northwest in the 
Northern and southwest in the Southern hemispheres, 
(ii) eastward-moving coherent dipolar structures of the 


type of equatorial modons, which do not appear in the 
“dry” adjustment. This experiment was repeated in a 
full sphere by Zhao et al. (2021) and they confirmed the 
original finding. The same adjustment scenario with baro- 
clinic structure led to generation of an MJO-like structure 
(Rostami and Zeitlin, 2020a). What was missed in all the 
aforementioned studies was the role of gradients of poten- 
tial temperature over the tropical zone. Thus, the two 
categories of equatorial adjustment that are presented in 
this study encompass simulations of the evolution of buoy- 
ancy over the equatorial zone. These experiments are (1) 
equatorial adjustment due to a large-scale localized nega- 
tive thickness anomaly, and (2) equatorial adjustment due 
to a large-scale localized positive buoyancy anomaly in 
the lower troposphere. Results show that both approaches 
produce almost similar results for large-scale eastward 
propagation that resembles the MJO-like structure. 


3.1.1 | Early stages of the adjustment 

Figure 3 presents the adjustment of a symmetric negative 
thickness anomaly that somehow represents a pressure 
anomaly in “dry” and moist-convective environments. 
Asymmetric perturbation would yield the emergence of 
Yanai waves owing to the response of the system to the 
initial depression. At the onset of thermo-geostrophic 
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adjustment, a quasi-isotropic inflow moves toward 
the equatorial center of the anomaly. For a large-scale 
anomaly, we also observe wind intrusion from the extra- 
tropical regions. Emission of a packet of fast and short 
inertia-gravity waves occurs in the initial stages. Then the 
Coriolis force causes a cyclonic circulation of the veloc- 
ity on both sides of the equatorial belt. In the presence 
of moist convection, these cyclones become pronounced 
and are converted into a quasi-equatorial modon. If the 
moist convection is weak, then these cyclones will propa- 
gate westward as Rossby waves. The other key difference 
between a “dry” and moist-convective environment is that 
detachment of long-wave BKW happens much sooner 
in a “dry” case. The response to adjustment in the “dry” 
case is similar to Gill’s mechanism (Gill, 1980); that is, 
Rossby-wave emission in the western sector and KW 
emission in the eastern sector of the localized anomaly. 
However, Figure 3 also indicates that a tipping point 
exists, where an eastward-propagating dipolar structure 
stems out in the presence of moist-convective diabatic 
processes as soon as the pressure anomaly reaches this 
threshold level (AH, /H ~ 0.07). For the case of a similar 
eastward-propagating equatorial modon in a barotropic 
structure, we need a greater pressure anomaly as the 
threshold value (Rostami and Zeitlin, 2019b). Let us recall 
that latent heat release leads to the intensification of a 
self-sustained negative pressure anomaly. The main tran- 
sition from Rossby wave to dipolar or quasi-equatorial 
modon occurs in the initial stages of adjustment. In 
baroclinic configurations, this transition occurs under a 
lower amplitude of pressure anomaly in comparison with 
barotropic structures (Rostami and Zeitlin, 2019a). Never- 
theless, fusion of the dipole with CCBKW and creation of 
a quadruple structure (Figure 4a,b) are purely a baroclinic 
effect. The condensation slows down CCBKW, and at the 
same time enhances the dipolar gyres in such a way that 
it can propagate eastward like an equatorial modon. The 
phase velocity of the CCBKW still has superiority over 
that of the vortex pair. Once the dipolar structure is gen- 
erated, it can continue its eastward propagation even in a 
“dry” environment, although dissipation and baroclinicity 
can slow down this propagation. 

As known, a negative pressure anomaly can generally 
be translated to a positive potential temperature anomaly, 
which is denoted as bı in the mcTRSW equations. We 
impose a similar large-scale localized disturbance with 
positive sign for buoyancy in the lower layer. Most afore- 
mentioned mechanisms can also be observed in the 
adjustment of potential temperature anomaly (Figure 5). 
Genesis of “hybrid structure” is associated with eastward 
advection of the buoyancy anomaly that is accompanied 
by the dipolar gyre circulation in a moist-convective envi- 
ronment. Figure 5 presents the evolution of this buoyancy 
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anomaly, which is a function of potential temperature, 
and the corresponding vorticity and condensation fields. 
Dipolar warm cores, which are surrounded by westerly 
and easterly winds, are prominent sources of the anoma- 
lous ascending motions that promote deeper convection. 
The emerging structure is similar to an equatorial modon. 
Such a similarity is evident in Figure A1, which presents 
the streamlines and velocity field of an asymptotic modon 
in stationary and comoving frames. Warm-front condensa- 
tion travels ahead of the eastward-propagating buoyancy 
anomaly (Figure 2). Formation of a “hybrid structure” 
can be seen in the initial stages of vorticity evolution 
(Figure 4a-c). Moist convection evidently plays a pivotal 
role at the center of this fusion in order to couple the dipole 
and the CCBKW (Figure 4c). Enhancement of the moist 
convection, which can be achieved by lower values of y, 
yields stronger attachment of the dipolar structure and the 
CCBKW;; otherwise there will be a weak westerly inflow 
in the gap between these two structures. Moist convec- 
tion intensifies these vorticity anomalies. Although at the 
onset of equatorial adjustment signs of vorticities are simi- 
lar, soon a weaker vortex pair with opposite signs arises in 
the upper troposphere and the whole system propagates 
slowly eastward. A weaker oppositely signed vorticity 
zone builds up on the eastern edge of the dipole due to 
the anticyclonic velocity field of the CCBKW. The quadru- 
ple is therefore zonally asymmetric. The two layers have 
opposite signs of vorticity, particularly outside the hybrid’s 
center in the eastern sector. In the western sector of the 
hybrid’s center, a “quasi-barotropic” dipolar structure is 
dominant. 

The intensity of condensation enhances abruptly at the 
very beginning of the equatorial adjustment (Figure 6), 
when condensation reaches its peak in the convergent 
zone. Then the spatial span of condensation shrinks 
with time while its intensity remains about the same 
(Figure 4c). The extreme condensation is accompanied 
by emission of inertia-gravity waves (IGW), which is 
presented by wave activity in Figure 6. This stimulated 
state, which is accompanied by dipole genesis, is not 
repeated during the whole evolution of excited struc- 
tures. A similar external forcing like intense large-scale 
localized heating over the warm pool in the lower tro- 
posphere, which is able to produce a similar velocity 
field, in principle, can create similar evolution. However, 
the existence of such a force in the atmosphere seems 
rare. Our key finding is that baroclinic adjustment in the 
lower troposphere in the presence of strong enough moist 
convection can create such a forcing in a self-sustained 
and self-propelled manner. The emerged dipolar structure 
(Figure 9), as a hybrid dipole-CCBKW, has a similar back- 
bone structure to that of the adjusted modon (Rostami and 
Zeitlin, 2019a). 
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Field of vorticity anomaly (colors) in the upper (a) and lower (b) layers and corresponding velocity field (arrows) at the 


initial stages of equatorial adjustment due to an initial negative thickness anomaly in the lower layer in a moist-convective environment. (c) 


Corresponding evolution of condensation (color) and baroclinic velocity field (arrows) in the lower layer. Time = 1, 4, 13, 24 [day], 


consecutively from left to right [Colour figure can be viewed at wileyonlinelibrary.com] 


3.1.2 | Late stages of the adjustment 

Initial and late stages of adjustment exhibit an eastward 
propagation (Figures 7-9). Eastward propagation of the 
buoyancy field, b4, also exists in a one-layer barotropic con- 
figuration, as seen by Kurganov et al. (2020). The buoyancy 
field acquires an off-equatorial advection at late stages. 
After the initial stages of the adjustment, we observe an 
expansion of the CCBKW (Figures 7 and 8), which demon- 
strates the main source of easterlies in the lower layer. 
Eventual detachment of the CCBKW leads to weakening 


of eastward-propagating westerlies at late stages (Figures 7 
and 8). Indeed, from Hovmöller diagrams we can infer 
that detachment of the CCBKW was never complete, with 
a significant easterly flow remaining in the frontal sec- 
tor of the dipolar westerlies, producing a highly conver- 
gent zone. It has already been proven theoretically and 
numerically that baroclinic modon solutions, as a “coher- 
ent” structure, do not exist in the 2RSW model (Rostami 
and Zeitlin, 2020b). The equatorial modon in this context 
can therefore be regarded as quasibarotropic, with a ten- 
dency for barotropization. This feature is compatible with 
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FIGURE 5 Evolution of different fields at the initial stages of equatorial adjustment due to an initial positive buoyancy anomaly in the 
lower layer in (a) “dry” and (b, c, d) moist-convective environments. Evolution of buoyancy (colors) and corresponding velocity field (arrows) 
in the lower layer in (a) “dry” and (b) moist-convective environments. (c) Corresponding snapshots of vorticity (colors) and velocity field 
(arrows). (d) Corresponding snapshots of condensation (colors) and baroclinic velocity field (arrows). Time = 1, 4, 13, 24 [day], consecutively 
from left to right [Colour figure can be viewed at wileyonlinelibrary.com] 
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the MJO that eventually ceases its eastward propagation. 
Both barotropic KWs and CCBKWs are interacting with 
the main dipolar structure in such a way that the former 
has about 10 interactions while the latter has at least one 
(Figures 7 and 8). Our results indicate that these interac- 
tions lead to loss of coherency and consequently decelera- 
tion of eastward propagation. The other factor that causes 
deceleration of eastward propagation is thermal relax- 
ation, which is introduced as Newtonian cooling in the 
model. Parameter y* controls the strength of this forcing 
in such a way that greater (smaller) values of y* decelerate 
(ameliorate; see Figure 1) the eastward propagation. 
During the elongation of the CCBKW, a secondary 
eastward-propagating dipole emerges from its western sec- 
tor (Figures 7 and 9). The first dipolar structure that lags 
behind eventually turns into westward Rossby waves. We 
observed a similar phenomenon in another experiment 
initialized with a positive buoyancy anomaly (Figure 9). 
While the whole hybrid system is moving eastward, a 
Rossby-wave tail has a tendency to move westward, which 


can cause eventual condensation, westward heat transfer, 
and interaction with circumnavigating BKWs (Figures 8 
and 9). The center of the “hybrid structure”, as the most 
convergent zone, exhibits the most condensation. An 
east-west asymmetry exists not only in condensation pat- 
terns but also in buoyancy evolution (Figure 9). 

One other remarkable finding of this study is that 
interaction of the circumnavigated BKW with the remain- 
ing eastward-propagating front at late stages leads to 
large-scale wave-breaking of the eastward-propagating 
front, as seen in the vorticity field (Figure 9) or Hovmöller 
diagrams (Figure 8). Full collision of this wave triggers 
westward-propagating Rossby waves in both main experi- 
ments (Figures 7 and 8). 

In the present article, we address the question of recur- 
rent generation of MJO-like structures. When BKW, after 
circumnavigating all around the Equator, interacts with 
a new large-scale positive buoyancy anomaly, bı, then 
this interaction triggers recurrent generation, with a new 
eastward-propagating “hybrid structure” (Figure 10). The 
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warm water in the western Indian Ocean, next to the 
East African Coastal Current (EACC), could be a plausible 
candidate for such a buoyancy anomaly. In this dual inter- 
action case, the critical amplitude of bı can be much 


smaller than that of the adjustment scenario with just 
the buoyancy anomaly alone (cf. left panel of Figure 5b). 
In other words, such interaction facilitates the triggering 
state of the new MJO-like structure. When the “hybrid 
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anomaly in a moist-convective environment. All right panels relate to eastward propagation of the MJO-like structure due to equatorial 
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structure” emerges from the thermo-geostrophic adjust- 
ment, both CCBKW and dipolar structure have intercon- 
nected eastward propagation; however, after about two 
weeks, CCBKW, due to its higher phase speed, produces a 
distinctive double rainband structure. The more extended 
the initial distribution of the bı anomaly, the later the 
occurrence of the double rainband structure. By approach- 
ing the suppressed stage of eastward propagation, the 
dipolar structure gradually loses its coherent structure 
and vorticity strength. After detachment, the BKW is 


accelerated due to less moist convection and transition to 
a (semi-)“dry” environment. The average phase speed of 
BKW in the Western Hemisphere is greater than that of the 
eastern one. 

The other result of this study is that, in addition to 
amplitude, the zonal scale of the buoyancy perturbation 
needs to be of approximately the same order as that of 
BKW to induce MJO-like eastward propagation. If the per- 
turbation scale is much smaller than that of BKW, then 
a fully adjusted, well-intensified vortex-pair structure 
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stages of equatorial adjustment in a moist-convective environment on a synoptic scale, and (b) large-scale meridionally elongated buoyancy 


perturbation. Time = 0, 10, 30, 50 [day], consecutively from left to right [Colour figure can be viewed at wileyonlinelibrary.com] 


cannot develop fully, propagate eastward, and couple with 
BKW (Figure 11a). What we see is like what Gill’s mech- 
anism predicts. Consequently, BKW acquires its normal 
phase speed and is not slowed down. The late-stage evolu- 
tion of a synoptic-scale buoyancy perturbation illustrates 
an eastward curved front. Evolution ofa large-scale merid- 
ionally elongated buoyancy perturbation in the lower 
layer (bj) illustrates that, if the meridional length-scale is 
much greater than that of BKW, then the modon loses its 
coherency. Consequently, the northern and southern vor- 
tices of the dipolar structure move to the northwest and 
southwest, respectively, and we see formation of a new 
dipolar structure on a smaller scale (Figure 11b). After the 
dissipation of the dipolar structure, b; reaches a ceased 
state of eastward propagation. 

We repeated a similar equatorial adjustment of a pos- 
itive buoyancy anomaly with realistic bottom topography 
(not shown) and found asymmetric behavior of the east- 
ward propagation due to the interaction of large-scale 
motions with land. During the initial days of adjustment, 
there is a meridional shift of buoyancy to the South- 
ern Hemisphere. This tendency can be ascribed to more 
enhanced moist convection in the southern sector of 
the Indian Ocean in comparison with the northern sec- 
tor, which is covered by land at latitudes above 15°. The 


dipolar component of the “hybrid structure” therefore 
rotates clockwise. As time goes on, buoyancy anomalies 
are dissipated and advected toward higher latitudes. A 
comprehensive study of adjustment scenarios, including 
land and topography, awaits further study. 


4 | COMPARISON OF THE 
RESULTS WITH THE MJO 


We summarize below some key features of the MJO 
episodes and their counterparts in our experiments. 


* Convective activity at the onset of MJO initiation 

A distinguishing feature of the MJO’s genesis 
in comparison with other stages of its evolution is 
an abrupt transition of convective activity from the 
intertropical convergence zone (ITCZ) to the Equa- 
tor (Chen et al., 2016). The strong dry air intrusion 
from the extratropical regions reaches its maximum 
equatorward wind speed during the MJO’s life cycle 
(Chen et al., 2016). This phenomenon is consistent 
with what we observe in the genesis mechanism 
based on the aforementioned adjustment scenarios 
(cf. Figure 3). 
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According to in situ and satellite observations, tran- 
sition from preconvective initiation to initiation stages 
of the MJO takes 2-20 days (Riley et al., 2011; Xu and 
Rutledge, 2015). Within this broad time span, daytime 
growth of the areal coverage and depth of shallow, 
nonprecipitating cumulus clouds exists (Ruppert and 
Johnson, 2015). Some observations illustrate precur- 
sor signals associated with MJO convection initiation 
about 3-10 days prior to the initiation date as anoma- 
lous ascending motion, heat fluxes (cf. Li et al., 2015; 
Xu and Rutledge, 2015), or enhanced humidity (Powell 
and Houze, 2014) in the lower troposphere. A remark- 
able coupling exists between moisture and convection, 
as seen from the wavenumbers and frequencies of 
the MJO episodes (Yasunaga and Mapes, 2012). This 
sort of preconvective signal in our simulations can be 
seen consistently from initial adjustment that is asso- 
ciated with a gradual increase of humidity or a tran- 
sition to deep convection, shown as abrupt growth 
of condensation within the initial days of adjustment 
(Figures 6 and 8). Enhanced convective activity is a nec- 
essary condition for evolution of MJO-like structures 
by equatorial adjustment, as seen in the comparison 
of “dry” versus moist-convective evolution (Figures 4c 
and 5) that occurs in the vicinity of the center of the 
depressed zone or positive buoyancy anomaly. Tran- 
sition to the onset of MJO initiation over the Indian 
Ocean is associated with deep cloudiness, enhanced 
convective activity, a series of quasistationary deep 
heating anomalies, and low-level moisture convergence 
(e.g. Hsu et al., 1990; Matthews and Kiladis, 1999; 
Kemball-Cook and Weare, 2001; Hsu and Lee, 2005; Seo 
and Kumar, 2008; Zhang and Ling, 2012; Straub, 2013). 
Figure 6 presents a peak level of condensation in the ini- 
tial stages of the phase-speed transition. A wacky look 
at this figure shows that two peaks exist in the gen- 
esis stages of the MJO-like structure: (1) at the very 
beginning of the adjustment, before formation of dipo- 
lar structure due to the high convergent zone, and 
(2) when the tower of elliptical condensation collapses 
toward the center and moves to the east, as shown 
in the first and second panels of Figures 4 and 5d. A 
repetitive pattern of a deep tropopause-penetrating con- 
vective process, which happens about one week before 
the main MJO peak in convection, has been observed 
via Cloud-Aerosol Lidar and Infrared Pathfinder Satellite 
Observations (CALIPSO: Genio et al., 2012). Our sim- 
ulations reveal a similar pattern, as seen, for example, 
in the Hovmöller diagram of condensation (Figure 8). 
Some other observations also show a similar rapid loft- 
ing of moisture into the middle troposphere at the onset 
of deep convection (Kiladis et al., 2005). It is worth not- 
ing that progression from shallow to deep convection on 


many time-scales or moisture transition is not unique 
to the MJO. It exists in almost all tropical convective 
systems (Mapes et al., 2006; Kiladis et al., 2009). 


Quadrupolar structure and enhanced dipolar structure 
Quadruple structure with an opposite-sign vorticity 
anomaly in the upper layer is known as one of the key 
observational facts of MJO events (Rui and Wang, 1990; 
Kiladis et al., 2005; Wang et al., 2019). The dominance 
of the rear dipolar structure is well documented in the 
literature (e.g. Lau and Phillips, 1986; Matthews, 2000; 
Zhang and Ling, 2012; Wang et al., 2019). Our results 
illustrate that at the onset of adjustment we observe a 
similar structure and in the early stages an opposite sign 
of vorticity in the upper troposphere. This structure is 
zonally asymmetric, because the western sector of this 
quadruple is occupied by the dipole, while the eastern 
sector is inhabited by the CCBKW (Figure 4). Conse- 
quently, the zonal velocity field toward the center of 
the “hybrid structure” is also opposite in the two lay- 
ers, particularly around the CCBKW. In observations, 
MJO events are accompanied by a forced component of 
hybrid Rossby-Kelvin type (cf. Hendon and Salby, 1994; 
Dunkerton and Crum, 1995). We see smooth detach- 
ment of the CCBKW in simulations of this study (cf. 
Figure 7). A similar pattern has been observed and 
analyzed by Sobel and Kim (2012). Our simulations 
reveal the nonuniformity of the “hybrid structure” and 
its deformation from genesis to suppression stages. One 
of these deformations is a zonal spread of the velocity 
field due to westward and eastward pulling of Rossby 
waves and KWs (Figures 7 and 9), which resembles the 
observed evolution of the MJO (Kiladis et al., 2005). 
All the simulations (cf. Figures 3 and 5b-d) illustrate 
that CCBKW is dominant in the east of the central con- 
vective envelope, which agrees with data analysis of 
tropical wave activity (Straub and Kiladis, 2003). 


Propagation speed versus Kelvin-wave phase velocity 
The slow eastward propagation speed of the MJO 
(x 5m-s71) may be its the most dominant feature that 
does not exist in other equatorial waves (cf. Weick- 
mann et al., 1985; Knutson et al., 1986; Wheeler and 
Kiladis, 1999; Wheeler and Hendon, 2004). As revealed 
by Figures 7 and 8, we observe two distinguished east- 
ward velocity fronts in the Hovmöller diagram of baro- 
clinic zonal velocity. The slower one is for the dipo- 
lar structure, which has the same order as the MJO 
phase speed, while the other indicates the CCBKW with 
phase speed ~ 12 m-s~!, which agrees with observations 
exhibiting a range of% 10-17 m-s! (e.g. Hendon and 
Salby, 1994; Wheeler and Kiladis, 1999; Straub and 
Kiladis, 2002). Initial conditions, such as the spatial dis- 
tribution of anomalies, thickness of layers, strength of 
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Newtonian cooling, and stratification of layers, 62/01, 
are factors that lead to slight variation of CCBKW 
phase speed. Intensive moist convection and latent 
heat release decrease the phase velocity of CCBKWs 
(Figure 3) at the onset of baroclinic adjustment, which 
corresponds to the active phase of MJO initiation. 
Slowed down CCBKW is one of the key components of 
“hybrid structure”. Observations also present that phase 
speeds of associated CCBKWs during the active phase 
over the Indian Ocean are smaller than those of sup- 
pressed phases over the eastern Pacific Ocean, where 
the MJO convection has almost ceased (Madden and 
Julian, 1972; Hendon and Salby, 1994; Milliff and Mad- 
den, 1996; Weickmann et al., 1997; Matthews, 2000; 
Sperber, 2003; Roundy, 2008). 


Precipitation patterns and westerly and easterly inflows 


The center of the “hybrid structure”, between the 
dipole and CCBKW, is the most convergent zone and 
is denoted as the main condensation zone, particularly 
in the initial stages (cf. patterns of condensation in 
Figures 4c, 5d, and 9), which is consistent with precip- 
itation in the MJO structure (Zhang, 2005). Regions of 
dense and depleted humidity have been observed to the 
east and west of the hybrid’s center (e.g. Sperber, 2003; 
Kiladis et al., 2005; Adamec and Wallace, 2015). Such a 
distribution of water vapor is consistent with the con- 
densation and divergence fields of the “hybrid structure” 
(Figure 4c). The typical swallowtail pattern of precipita- 
tion, which is common for the MJO during active phases 
(Zhang and Ling, 2012; Adamec and Wallace, 2015), is 
consistent with snapshots of condensation (Equation 9). 

Low level (~850-hPa) westerlies and easterlies in 
the western and eastern sectors of the “hybrid struc- 
ture” are frequently observed characteristics of the MJO 
(Madden, 1986; Sperber, 2003; Zhang, 2005). The MJO 
includes a stiff wind over an area roughly 1,500 km 
in latitude and 4,500 km in longitude. The MJO has a 
great impact on both tropical moist-air circulation and 
midlatitude intraseasonal environment variables (Vitart 
and Molteni, 2010). Figures 4,7, and 9 indicate a sim- 
ilar pattern of westerlies and easterlies with the same 
order of magnitude. Our results reveal that, by dimin- 
ishing the strength of moist convection, there will be 
a gap between the dipolar structure and CCBKW that 
is filled up by weak westerlies (Figures 7 and 9). In 
this case the detachment is accompanied by a distinc- 
tive double rainband structure. One of these structures 
moves with the CCBKW, so its phase speed is in the 
range of that of CCBKWs. Some similar precipitation 
patterns have been documented by Chen et al. (2016) 
and Zhu and Li (2017). Let us note that condensation in 
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this version of the mcTRSW model is synonymous with 
precipitation and presents large-scale deep convection. 


Potential temperature anomaly in MJO episodes 


A steadily daily-averaged increase in sea-surface 
temperature (SST) prior to the active phase of the MJO 
and gradual slowing down of the SST during eastward 
propagation have been reported in observational stud- 
ies (cf. Vialard et al., 2008; Matthews et al., 2014; Moum 
et al., 2014; Chen et al., 2016). SST is not included in our 
stand-alone version of the model directly. If we assume 
that the buoyancy field in the initial step is affected by 
the SST, then Figures 5 and 9 will illustrate the evolution 
of the localized disturbance of atmospheric buoyancy 
in the lower layer. Figure 9 apparently shows that there 
is a lag between eastward propagation of the convec- 
tive zone, illustrated as condensation, and the potential 
temperature anomaly. A similar lag has been observed; 
for instance, Hendon and Glick (1997) reported that 
spatially coherent SST anomalies propagate eastward in 
the Indian Ocean along with the large-scale convective 
anomaly, but with a 1/4 cycle lag. 


Sub/mesoscale convective envelopes of high buoyancy at 
late stages 


One of the dominant processes that has a remark- 
able impact in destabilizing the eastward propaga- 
tion and genesis of cyclones at late stages is colli- 
sion of circumnavigated CCBKW into it. As is clearly 
visible through Hovmöller diagrams (Figures 7 and 
8), this destabilization not only diminishes the total 
eastward-propagating time span but also engenders a 
few convective events with smaller spatial scale and 
longevity that propagate westward, as is easily vis- 
ible at late stages of Hovmöller diagrams. In real 
MJO events, although detachment of smaller-scale 
low-frequency convective envelopes has been observed 
from the onset of the MJO’s active phase as detachment 
of a westward-propagating cloud cluster of order 100 km 
(Leary and Houze, 1979) from a large-scale super clus- 
ter (Nakazawa, 1988), at late stages the number of 
envelopes is increasing and some of these small con- 
vective envelopes can merge and form tropical storms 
in the western Pacific (Straub and Kiladis, 2003). The 
formation of similar low-frequency tropical cyclones in 
the eastern Pacific is ascribed as more likely due to 
MJO-related episodes with lower-tropospheric westerly 
anomalies (Maloney and Hartmann, 2000). According 
to Maloney and Hartmann (2000), westerly (easterly) 
equatorial wind anomalies in the lower troposphere 
are accompanied by enhanced (suppressed) convection 
over the eastern Pacific hurricane region. 
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e Recurrent generation and episodic versus periodic char- 
acter of the MJO 

Each spatial distribution and amplitude of bı gen- 
erates a “hybrid structure” that has its own longevity, 
appearance, and periodicity of the MJO-like structure 
(Figure 11), even though the underlying mechanism for 
each episode is the same. The intensity and extension 
of diabatic processes associated with the coupled struc- 
tures in each episode and consequently the amount of 
deceleration of BKW are not exactly the same. Thus, it 
is no wonder that the MJO is not strictly an oscillation 
but an episodic severe weather event (Hartmann and 
Hendon, 2007). 

The aforementioned mechanism for the recurrent 
generation of MJO-like structure (Section 3) also gives a 
plausible explanation for various studies that have doc- 
umented the connection of BKW to the next MJO event 
and its role in activating convection of the next cycle of 
the MJO once BKW enters over the warm water in the 
western Indian Ocean next to the EACC (cf. Knutson 
and Weickmann, 1987; Kikuchi and Takayabu, 2003). 


5 | CONCLUSION AND 
DISCUSSION 


By means of a new pseudo-spectral multilayer mcTRSW 
model in a full sphere that uses spin-weighted spheri- 
cal harmonics, we extended the MJO mechanism pre- 
viously proposed by Rostami and Zeitlin (2020a). In 
this study, potential temperature variation, which was 
the missing part of the mcRSW model used by Ros- 
tami and Zeitlin (2020a), is implemented in the mcTRSW 
model in a full sphere. In contrast to the classical 
setup of shallow-water models, gradients of the verti- 
cally averaged potential temperature and density are vari- 
ant in the mcTRSW model. Therefore, in addition to 
studying the baroclinic evolution of the relaxation of 
large-scale localized warming due to a negative pres- 
sure anomaly over the equatorial zone, we could exam- 
ine the explicit evolution of a buoyancy anomaly. We 
found an equatorial adjustment beyond Gill’s mecha- 
nism (Gill, 1980) for the genesis and backbone structure 
of an eastward-propagating MJO-like structure. Accord- 
ing to this theory, an eastward-moving structure can be 
generated in a self-sustained and self-propelled manner 
whenever the spatial scale and amplitude of the buoy- 
ancy anomaly or negative pressure anomaly or both reach 
a critical threshold state in the lower troposphere in the 
presence of moist convection. The excited system has a 
“hybrid structure” that is constituted by a combination of 
a “quasi-equatorial modon” and a CCBKW. The latter has 
greater phase speed than the former. This coupling lasts for 


an intraseasonal period. This mechanism captures most of 
the crudest features of the MJO during the genesis, early, 
and late stages of its evolution in the lower and upper 
troposphere, as seen in the convective activity, quadrupo- 
lar structure, vorticity fields, wave activity, phase speed, 
westerly and easterly inflows, condensation pattern, and 
intraseasonal lifespan in an obvious manner. Inclusion of 
the buoyancy term in the mcRSW model enhances the cou- 
pling between dynamics and convection and hence yields 
a more realistic evolution of the “hybrid structure”. We 
examined the same adjustment scenarios with real bot- 
tom topography as well, and we found that the crudest 
features of the aforementioned adjustment mechanism for 
generation of the “hybrid structure” are robust. 

Another important question that we addressed in the 
present article is a possible mechanism for the recurrent 
generation of the next cycle of the MJO. At the late stages 
of evolution of the “hybrid structure”, the BKW becomes 
less strongly attached to dipolar structure and convec- 
tion. Thus, it finally becomes detached and continues 
its circumnavigation all around the Equator with greater 
phase speed in (semi-)dry air. When the circumnavi- 
gated BKW, with intraseasonal periodicity, interacts with 
a new disturbance of the large-scale buoyancy anomaly, 
a new eastward-propagating “hybrid structure” emerges 
from thermo-geostrophic adjustment. The warm water in 
the western Indian Ocean, next to the EACC, can provide 
such large-scale buoyancy anomalies. Furthermore, the 
higher phase speed of the CCBKW with respect to that of 
the dipolar structure may yield double condensation fronts 
during the separation stage in the presence of moist air. 

In the present study, we also found a mechanism to 
explain the small-scale turbulence that appears in the late 
stages of simulations; that is, collision of the circumnav- 
igating convectively coupled BKW with a slow-evolving 
remaining part of the perturbation. This collision is accom- 
panied by subsequent destabilization of the vorticity shear 
and/or convective small-scale instabilities of an elon- 
gated potential temperature front. It is worth mention- 
ing that collision of the circumnavigated BKW with a 
large-scale positive buoyancy anomaly can trigger new 
eastward propagation of the “hybrid structure” that may 
explain the intraseasonal periodicity of MJO episodes. Suc- 
cessive interactions of the eastward-moving “hybrid struc- 
ture” with IGW, barotropic KW, and also thermal cooling 
have a destructive impact on coherency and a decelera- 
tion effect on the phase speed of the “hybrid structure” 
that lead to the generation of many small-scale distur- 
bances. The aforementioned adjustment scenarios could 
also explain the appearance of the double rainband struc- 
tures that have been observed in some MJO episodes as 
detachment of CCBKW from the dipolar structure in a 
weak moist-convective environment. 
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The observed mechanism in our model could provide 
some supplementary explanations for some existing the- 
ories. In this respect, it encompasses some theories that 
previously seemed divergent. For instance, the “hybrid 
structure” supports the existence of an equatorial modon 
or planetary-scale solitary Rossby-vortex pair in the MJO 
(cf. Yano and Tribbia, 2017) as a large-scale tropical cir- 
culation that is nondivergent to leading order. The dipolar 
structure can maintain its coherency and eastward prop- 
agation in the absence of dissipation and baroclinicity, 
even in a “dry” environment. Moist convection is a nec- 
essary condition for the genesis of the “hybrid structure”, 
at least in the early stages of equatorial adjustment. Moist 
convection in this mechanism yields slowed-down zonal 
propagation of CCBKW as well. By illustrating that a com- 
bined “hybrid structure” yields similar slowed-down prop- 
agation of CCBKW via moist convection, the results of the 
present study support theories that include a slowed-down 
Kelvin wave as a core component of the MJO’s structure, 
e.g. Kim and Zhang (2021), who produced a slowed-down 
Kelvin wave by including momentum damping in the 
linear equatorial shallow-water equations in a “dry” envi- 
ronment. Although the results in this article reveal high 
degrees of convergence with theories that include a mixed 
Kelvin-Rossby wave structure (e.g., Wang et al., 2016 
and references therein), there is a hermeneutic difference 
between the quasi-equatorial modon with eastward propa- 
gation and Rossby waves with westward propagation in the 
western sector of the MJO-like structure. The vortex-pair 
structure in our model emerges as a modon-like struc- 
ture beyond Gill’s mechanism when a buoyancy anomaly, 
depressed pressure, or a combination of both reaches 
a triggering state, while the theories based on mixed 
Kelvin-Rossby wave structure follow Gill’s mechanism. 
The vorticity of the vortex pair in a modon-like structure 
is more intensified than that of a Rossby wave in Gill’s 
mechanism. Capturing the full nonlinear evolution of the 
variables is another difference of our proposed mecha- 
nism from those of linear wave theories (cf. Fuchs and 
Raymond, 2017) that helps us to simulate more features 
of the MJO-like structure, like a well-intensified vortex 
pair, for the reason that, according to basic knowledge of 
wave theory, the phase speed of an equatorial modon, in 
its ideal form, is outside the range of speeds for a lin- 
ear Rossby wave of any zonal scale (cf. Osborne, 1991; 
Grimshaw and Iooss, 2003). Namely, Rossby waves need 
to be identified in the wavenumber-frequency phase space 
outside the range of linear wave dispersion to propagate 
eastward (Yano and Tribbia, 2017). Moreover, in the Trop- 
ics, the magnitude of the Rossby number of atmospheric 
motion is typically one order greater than in mid and 
high latitudes; therefore, the nonlinear effect of advec- 
tion becomes more important (Maicun, 1987). In moist 


Royal Meteorological Society 


wave theories (cf. Adames and Kim, 2016) eastward prop- 
agation is due mainly to meridional advection of mean 
moisture by MJO winds. In such moist wave theories, a 
prognostic moisture anomaly exists in the parameteriza- 
tion of the boundary layer and a passive Gill-like response. 
However, in the proposed model in this study there is 
no need for prognostic moisture. Moisture-fed convection 
provides a mechanism for the “hybrid structure” to be 
excited and maintained. This mechanism differs funda- 
mentally from that of moisture-mode theories, in which 
no similar “dry” dynamical basic structure exists. The 
barotropic equatorial modon and BKWs have innate east- 
ward propagation even in a “dry” environment. By grow- 
ing baroclinicity in the equatorial modon, we observe loss 
of coherency and an arrest of eastward propagation (Ros- 
tami and Zeitlin, 2020b). 

This study aims at understanding the fundamental 
dynamics of the MJO. Thus, some features of the MJO, 
such as details of the vertical structure, full boundary 
effects of the maritime continent, effects of bottom topog- 
raphy, direct effects of SST, and influence of radiative 
transfer have not been captured by this hierarchy of the 
model; nevertheless, these features can be investigated by 
more complex versions of the model. This model, with its 
intermediate complexity, is also the dynamical core of Aeo- 
lus2, one of the atmospheric models of the Potsdam Earth 
Model (POEM). This model has the capability to capture 
thermal and dynamical effects of orography and radiative 
transfer distinctively with physically reasonable asymp- 
totic limits. Via relaxation of the potential temperature of 
the air to SST, the model is able to be linked to SST vari- 
ation of an oceanic model. As model Aeolus2 is based on 
minimal parametrization, it can be used for testing general 
circulation models too. 
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APPENDIX A. SUMMARY OF ANALYTIC 
SOLUTION OF ASYMPTOTIC MODON 
STRUCTURE 


Following Rostami and Zeitlin (2019a), let us recall that 
the equation for the streamfunction of the equatorial 


Royal Meteorological Society 


modon is derived from RSW equations in the equatorial 
p-plane with no dissipation: 


aviv: W+fhykaAv+gvh=0, 
ðh +V- (wh) =0, 


(Al) 


where V = (0,, 0,), v = (u, v), u and v are zonal and merid- 
ional components of the velocity, h is geopotential height, 
p is the meridional gradient of the Coriolis parame- 
ter, and k is a unit vertical vector. A small pressure 
perturbation (A > 0) can be included in the equations 
using parameter A: h = H(1 + An), where H is unperturbed 
geopotential height. Spatial, velocity and time-scales are 
(x,y) ~ L, (u,v) ~ V, t~ L/V. Under the hypothesis that 
the characteristic velocity is much smaller than the phase 
velocity of the barotropic Kelvin waves, c = ygH, V << 
gH , and gHA/V? = @(1), hence, the nondimensional 
equations take the form 


aviv: W+BhykAv+Vn=0, (A2) 


Man +v- Vn) +(1+ dn) V-v=0, (A3) 


where p = pL? /V, and v = Vo + Av, +-+- . In the leading 
order in A, Equation A3 gives V - vo = 0, the motion is non- 
divergent, and uo = —dy;y, vo = Oxy. The streamfunction 
can be extracted by cross-differentiation of the zonal and 
meridional momentum equations: 


Vwi + Tw, Vy) + pyx = 0, (A4) 


where J is Jacobian. The modon solutions are obtained 
under the hypothesis of steady motion with constant zonal 
velocity U, by assuming a linear relationship between the 
absolute vorticity and streamfunction in the comoving 
frame, which gives inhomogeneous Helmholtz equation. 
It is solved by separation of variables in polar coordinates 
in terms of Bessel functions, first in the outer domain 
under the decay condition, and then in the inner domain, 
and matching the inner and outer solutions across a circle 
of given radius a in the plane. The solution is in the form: 


_ _Ua 
K, (pa) 


Wext = K,(pr)sin @, r> a, 
Ton [hov SLU uv?)| sind, r<a, 

(A5) 

where Jı and KĶ; are ordinary and modified Bessel func- 

tions of order one, respectively, p is real, and p? = B/ U, 

so U > 0, meaning that the motion is eastward. For each 

pair (a, p) there exists a series of eigenvalues v arising from 

matching conditions, of which the lowest corresponds to a 
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FIGURE Al Streamlines and velocity field of an asymptotic modon in stationary (left panel) and comoving (right panel) frames with 
U = 0.1. Dashed circle: separatrix of radius a = 0.5 [Colour figure can be viewed at wileyonlinelibrary.com] 


dipole, to be called the asymptotic modon. The next eigen- 
value gives a quadrupolar solution called the “shielded 
modon,” and so on. The higher eigenvalues correspond to 
multipolar vortices with significant velocity shears in the 


core, which are more sensitive to dissipation. Centering 
the solution in the equatorial #-plane leads to a cyclonic 
pair. The streamlines of the asymptotic modon according 
to Equations A5 are presented in Figure A1. 


ASUDII'T SUDO OANAID a[quortdde oy) Aq pouroAo’d ae SANIE YO ‘asn Jo SMI Joy Á1eIQIT IUUQ ÁM UO (SUOTTpUOD-pue-sULIO]/WOS ÁM AreAQHouT[Uo//:sdyYy) SUOTIPUOD pur SUAL IY? 29g *[ZZOZ/T1/OT] wo Arerqr] IUUQ ÁM ‘PULI oueIyD0_D Aq ggep !b/ZOOT OT /Lop/mos Kop Arerqrjaurfuo'sjouuy//:sdyy wory papropumod ‘0 ‘XOL8LL+1 


